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Abstract 

We study the vector-vector interaction within the framework of the hidden gauge 
formalism for the channels with quantum numbers Charm C = and Strangeness 
S = in the energy region around 4000 MeV. By looking for poles in the complex 
plane we find three resonances that could be identified by the mass, width and 
quantum numbers with the Y(3940), Z(3940) and X(4160), these poles appear with 
^j-' \ isospin 1 = and J PC = ++ , 2 ++ and 2 ++ respectively. Whereas the Y(3940) 

and Z(3940) are coupled more strongly to D*D*, the X(4160) is basically a D*D* 
molecular state. Another two extra resonances appear in our approach with / = 0, 1 
and J PC = l + -,2 ++ which are not found in the PDG with masses M = 3945,3912 
qs^ ' MeV and widths T = 0, 120 MeV respectively. 

O ' 

: 1 Introduction 

The charmonium spectroscopy has been recently pushed forward with the unexpected 
discovery of many new charmonium-like resonances in the B factories at SLAC, KEK 
and CESR. The B-factories, originally constructed to test matter-antimatter asymmetries 
or CP-violation, within or beyond the standard model, discovered a number of interesting 
charm and hidden charm mesons which do not seem to have a simple cc structure. The first 
of these XYZ states is the V(3872) that was observed by the Belle Collaboration as a narrow 
peak near 3872 MeV in the 7i + 7i^J/tp invariant mass distributions in B~ — > K~7t + tt~ J/ip 
decays [H [21 [31 H] . Belle also observed that the rate of the decay V(3872) — > 7r + 7r~7r° J/ip 
was comparable to 7r + 7r~ J ftp [5J, with the 3ir or 2ir coming from uj or p which would imply 
the V(3872) is a mixture of both 1 = and / = 1 |6] (see however [3 [9] for alternative 
explanations of that ratio). The close value of the V(3872) mass to the sum of the masses 
m D o + rriD*o = 3871.81 ± 0.36 MeV led to consider the V(3872) as a molecule-like bound 
state of a D° and a D*° meson, and much work on the properties of these possible systems 
has been done jSl C3 EB ED H21 E31 H3 • 
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The next XYZ states, that we will consider in this paper, are the X(3940), the F(3940), 
the X(4160) and the Z(3930). The X(3940) was observed in the double- charmonium 
production reaction e + e" — > J/ip + X with mass M = 3943 ± 8 MeV and width T < 52 
MeV [15J. After that Belle also observed a D*D* mass peak in the e + e~ — > J/tpD*D* 
reaction [16J. Whereas the X(3940) — > DD* has been observed, there is no signal for 
the DD or the ujJ/ip decays. Because of that and the fact that the t] c {1S) and r] c (2S) 
were also produced in double-charm production, it was believed that the X(3940) could 
have J PC = h , being a S^^Si charmonium state (rj"), but there are problems with this 
assignment since in this case the X(3940) should have a mass ~ 4050 MeV or even higher 
[T7] . Thus, it seems very unlikely that the X(3940) is a cc state. 

Belle has observed a state with a decay mode V(3940) — > uJ/ip in5^ KuJ/ip decays 
[T8] and Babar has confirmed it [19] , although the values for the mass and width reported 
by Babar are smaller than Belle's values (M = 3943 ± 17 MeV and T = 87 ± 34 MeV is 
reported by Belle and M = 3914.3±3;g MeV and T = 33±§ 2 MeV by Babar). In principle, 
the mass and width of the V (3940) suggest a radially excited P-wave charmonium state but 
then the x cl (2P) — > DD* would be the dominant decay mode and has not been observed. 
Hence, it seems very unlikely that the reported V(3940) and X(3940) represent different 
decay modes of the same state [20] . The J PC assignment in the case that the V(3940) were 
a charmonium state is not very clear. Indeed, for a charmonium state with J PC = h 
(r] c ) the mass is a little low and for J PC = ++ (x' c0 ) the mass is too high [20] . also if it 
had a simple cc structure one would expect that the open charm decay modes would be 
dominant, and the other, uJ/ip, negligible. Also, the large partial width V(3940) — > u.J /ip, 
estimated above 1 MeV [21], is quite larger than the measured partial widths for any of 
the observed hadronic transitions between charmonium states. 

Belle has recently observed a mass peak in the D*D* system in the e + e~ — > J/ipD*D* 
reaction [TB]. This state with a mass of (4156 ± 29) MeV and a width of T = 139^ 3 
MeV has been called A(4160). The production mechanism ensures that it has C — +. 
The known charmonium states seen from a e + e~ — > J/ipD*D* reaction have J = 0, thus, 
although this state could be identified with a 3 1 Sq (r/") or 4 1 5o (r)"') charmonium state, 
the mass predicted in the first case would be smaller (4050 MeV) and higher (4400 MeV) 
in the second case [17] . Recently, the CDF Collaboration at Fermilab has announced a 
narrow peak near the J/i/xf) threshold, which is designated as V(4140), observed in the 
B + -> J/ip(f)K + decay, it has a mass M = 4143 ± 2.9(stat) ± 1.2(syst) MeV and a width 
r = ll.7tH(stat) ± 3.7(syst) MeV [22]. The width observed is quite different from the 
one reported by Belle for the A (4160) [TB], suggesting that one is taking about a different 
state. Another point is the structure and quantum numbers. There are some predictions: 
In [23] the authors solved the Schrodinger equation from the potential obtained using 
effective Lagrangians and they found molecular solutions for V(3930) and V(4140) with 
J p = + ,2 + , concluding that the V(3930) and the V(4140) are molecular partners. In 
[24] , the authors assume that the V(3940) and the V(4140) are hadronic molecules with 
quantum numbers J PC = ++ or 2 ++ whose constituents are the charm vectors D*D* for 
the V (3940) and D* s + D* s ~ for the V (4140) and they calculate the decay rates of the observed 
modes V(3940) J/^juj and V(4140) J/ip for the case J PC = ++ . The coupling 
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constants are determined by means of the compositeness condition [24] and the results for 
these decay modes support the molecular interpretation of the Y(3930) and the K(4140). 
In [25l [26] the authors use QCD sum rules to evaluate the mass of a possible mesonic state 
that couples to a molecular D* + D*~ current. Whereas in [26] a mass M = 4.14 ±0.09 GeV 
is found, concluding that it is possible to describe the F(4140) as a D*D^ molecular state, 
in [25] they find a larger value of the mass M = 4.43 ±0.16 GeV. In [27] the authors suggest 
possible assignments for the charmonium-like "X, Y, Z" states, for example the X(4160) 
could be assigned to a Xco{3P) state. Also in [28] the authors find that the P-charmonium 
x"j(J — 0, 1) associated to the K(4140) is problematic because of the upper limit of the 
branching ratio of the decay Y(4140) — > J /ip<f) computed of the order of 10~ 4 — 10 -3 . In 
[29] a tetraquark structure is suggested for this state, while in [30] it is argued than the 
peak at this energy is just a <fiJ /ip threshold effect. 

The Z(3930) is also reported by Belle as a peak in the spectrum of DD mesons produced 
in 77 collisions, with mass and width M = 3929 ± 6 MeV and T = 29 ± 10 MeV. In the 
production process the two photons can only produce DD in a ++ or 2 ++ states. The 
Belle measurements favors the 2 ++ hypothesis, making the assignment of the Z(3930) to 
the 2 3 P 2 (Xc2) charmonium state possible [21], following the arguments given in |27j . 

In this paper we will propose a theoretical explanation on the nature of some of these 
XYZ states, providing structure and quantum numbers for them. Our work is based on 
the hidden gauge symmetry (HGS) formalism for the interaction of vector mesons, which 
was introduced by Bando-Kugo-Yamawaki [31], [321, [33] . This hidden gauge formalism has 
been applied in [34] for the pp interaction, giving rise to two bound pp states that could be 
identified with the / (1370) and / 2 (1270). After that, the radiative decay of these states 
into 77 was studied in [35], obtaining results in good agreement with the PDG [36]. The 
work of [34] was later extended to SU(3) [37], where several states were found that could 
also be identified with some of the PDG. Recently, the pD*, uD* interaction was studied 
within the same formalism [38], generalizing to SU(4) the basic hidden gauge Lagrangians, 
but breaking the SU(4) symmetry in the vector exchange diagrams. Some states were 
found which could be identified with the .02(2460) and the _D*(2640), giving a prediction 
of the quantum numbers for the last one I(J P ) = 1/2(1 + ) in the case the correspondence 
was correct, and providing a reasonable explanation for the small width of this state. Also, 
a new D Q state was predicted with J p = + , a mass close to 2600 MeV and width of about 
60 MeV. 

The present paper follows the steps of [38] , investigating possible vector - vector states 
of hidden charm, using a unitary approach in coupled channels including all the channels 
involved for Charm C = and Strangeness S = 0. 
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2 Formalism for VV interaction 



2.1 Lagrangian 

Our starting point is the Lagrangian, which involves the interaction of vector mesons 
amongst themselves, coming from the formalism of the hidden gauge symmetry (HGS) for 
vector mesons [3H [32], [39] 

Cm = ~\{V,uVn , (1) 
where the symbol () stands for the trace in the SU(4) space and is given by 

= dfjVy - d u V, - ig[V^ V v ] , (2) 

with g given by 

My 

9 = W , (3) 

and / = 93 MeV the pion decay constant. Using the value of g in Eq. (j3J) is one of the 
ways to account for the KSFR relation [40J which is tied to vector meson dominance [H] • 
The vector field is represented by the 577(4) matrix which is parametrized by 16 vector 
mesons including the 15-plet and singlet of S77(4), 



p- -£ 2 + % K*° D*~ 
K*~ K*° <p D*- 

V D*° D* + D* + J/^p J 



(4) 



where the ideal mixing has been taken for lu, <ft and J/ip. The interaction of Cjjj gives rise 
to a contact term coming from [V^, V^jfV^, V u ] 

4// = jiV^WV - V v V^V v ) , (5) 
depicted in Fig. d]a), and on the other hand it gives rise to a three vector vertex 

£ff = iS«^-Wn , (6) 

depicted in Fig. [T]b). This latter Lagrangian gives rise to a VV — > VV interaction by 
means of the exchange of one of the vectors, as shown in Figs. H]c), d). 

The SU (4) structure of the Lagrangian allows us to take into account all the channels 
within SU(4) which couple to certain quantum numbers. In the present work we shall 
present results for the case of all the channels which couple to D*D* or D*D*. The 
formalism is the same used in [31] and [38J. Some approximations were made there which 
make the formalism handy and reliable, by neglecting the three-momentum of the vector 
mesons with respect to their masses. It is interesting to see that with this approximation 
one obtains [3H] from the hidden gauge approach the chiral local Lagrangians which are used 
to study the interaction of pseudoscalar mesons among themselves and the pseudoscalar 
mesons with vector mesons and with baryons [121 H3], HH HH] • 
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Figure 1: Terms of the Cui Lagrangian: a) four vector contact term, Eq. b) three- 
vector interaction, Eq. (151); c) t and u channels from vector exchange; d) s channel for 
vector exchange. 

D* + (h)\ /D*+(k 3 ) 
D*-(k 2 )/ \ D*-(fa) 



Figure 2: Contact term of the D* + D* interaction. 
2.2 Four-vector contact interaction 

The channels that we are interested in are those with Charm C = and Strangeness 5 = 0, 
they are, in the case of I = 0: 

D*D*(4017), D*D*(4225), K*K*(1783), pp(1551), 0^(1565) 

00(2039), J/V»J/V>(6194), cjJ/^(3880), 0J/^(4116), ^0(1802) 

where the magnitude between parenthesis is the sum of the masses of the two meson 
involved, and for 1 = 1: 

D*D*(4017), K*K*(1783), pp(1551), pw(1558), pJ^(3872), p0(1795). 

We are not interested in the case of / = 2 which was considered in [37J and where no bound 
states or resonances were found. 

Consider now the D* + D*~ — > D* + D*~ reaction, see Fig.[2j In order to get the amplitude 
we use the Lagrangian of Eq. (j5J) obtaining 

C +D ^o* +D *- = 2^^^" + ef^^"- 2e^e^n , (7) 

where the indices 1, 2, 3 and 4 correspond to the particles with momenta ki, &2, &3 and k± 
in Fig. [2j It is straightforward to write down all amplitudes for the other channels. 

In the approximation of neglecting the three-momenta of the vector mesons, only the 
spatial components of the polarization vectors are nonvanishing, the orbital angular mo- 
mentum is L = 0, and then one can obtain easily spin projection operators [34J into spin 
0, 1, 2 states, which are given below: 
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P (0) = \e,eV 

= ^(e^eV-e^eV) 

P (2) = {^( V /6 v + V /^)-^/6/}. (8) 

where the order 1, 2, 3, 4 in the polarization vectors is understood. We can then write the 
combination of polarization vectors appearing in Eq. ([7]) in terms of the spin combinations 
and thus we obtain the kernel of the interaction which will be later on used to solve the 
Bethe-Salpeter equation. However, it is practical to construct the isospin combinations 
before the spin projection is done. 

Recalling that we have isospin doublets (-D*°,D* + ) and (D*~,D*°), the / = and 1 
combinations are written as 



\D*D*I = 0,I 3 = 0) = -L\D* + D*-) + -^\D*°D*°), 
l > / ^ I /T /, 

!£>*£>*, J = 1, lg = 0) = -^=\D* + D*-) - ^|L>*°.D* ). (9) 

We then find the amplitudes in the isospin base by forming linear combinations of the 
amplitudes in the particle base weighted by the Clebsh-Gordan coefficients as given in 
Eq. ®, 

f ('=°> = 3 g 2 (e^e v e v + e^e^ - 2e^e»e») , 

t {1=1) = g 2 {t^e» + t^e v -2t^e u e») . (10) 

These amplitudes, after projection into the spin channels, give rise to the following kernels 
(potential) for / = 0, 

f( j=o,j=o) = +6g2) 

f (J=0,J=l) = +9g 2 ; 

t (/ =0 ,J=2) = _ 3 ^2 ; (n) 

and also for the case of I = 1, 

t (/=l,J=0) = 2 ^2 ? 

t (/=l,J=l) = 3g 2 ^ 

t (/=l,J=2) = _^2_ (12) 

In the same way, we proceed to calculate the contact term projected in isospin and also 
spin between two of anyone of the involved channels. 



6 



^* + (/ci,ei) 



Figure 3: three- vector vertex diagram. 



2.3 Vector meson exchange terms 

Continuing with the D* + D*~ — > D* + D*~ reaction, now we want to calculate the amplitude 
of the first diagram in Fig. HI In order to do that, we follow the same procedure as 
in [38], starting off the three- vector vertex which is given by the Lagrangian of Eq. ([6]) 
and is depicted in Fig. [3] in the case of D* + — ► D* + p°. The amplitude of Fig. [3] can be 
automatically calculated rewriting the Lagrangian of Eq. (jfJJ) as, 

= ig({V^d u V^-d u V^)V u ) . (13) 

As it was explained in [38], within the approximation of neglecting the three-momenta of 
external vectors, the V u corresponds to the exchanged vector, simplifying the calculation. 
As already mentioned, this corresponds to the consideration of only the s-wave. 
The vertex function corresponding to the diagram of Fig. [3] is given by 

= ^{h + hU lu elef (14) 

With this basic structure, and considering all the particles involved in the exchange, we 
can readily evaluate the amplitude of the first diagram of Fig. H] to obtain 



12 11 

t%V +D ,-^ D , +D ,- = -- 9 2 (j^ + Jp^ kl + ^ ' ^ + ^ ^^^4 ( 15 ) 

In the vector exchange diagrams the exchange of one heavy-vector meson becomes 
supressed with respect to the exchange of a light vector by the ratio M|/M^ with L = u>, p 
and H = J /ip. To project in s-wave one must do the following replacements: 

, , s-Mf- Ml 

h ■ k 2 = — 



2 

k\ ■ ks = kik% — p- q 



u,n (s + Ml - M|)(s + Ml - Mi 



4s 
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Figure 4: Vector exchange diagrams for D*D* — > D*D*. 

where — > means the projection over s-wave, and k\ = (kx,p), k 2 = (k®, —p), ^3 = (k^q), 
k\ = (&4, —q) and Mj, with i = 1,4, is the mass of each particle, and so on for the other 
products of momenta. Thus, Eq. (115p can be written as 



.(„) ff 2 (2M 2 M 2 + Mj^(M^ + M 2 ))(4M 2 » - 3 g )e lM e 2 ,e^ 



The isospin projections taking into account the three diagrams of Fig. H] and Eq. 
give us 

{ex j=o) 9 2 (2MZM} + M 2 ^(3M 2 + M;))(4M^, -3^e^ 



^=1) _ 9\2MjMj + Mj^(-M j + MgK_4 M^ - 3s)e M e 2 „e^ 



,(ex,l=l) & \ m—p 1 "~j/if>\ —a; 1 —p//\ -t^- f-|7\ 



Now using the equations for the spin projections we can split the terms into their spin 
parts and we obtain 

, (e ,,/=o,5=o,i, 2 ) 9 2 (2 M 2 M 2 + M^(3 M 2 + M 2 ))(4 M|,. - 3 s) 

, {e ,,/=i,5=o,i,2) ^ 2 (2 M 2 M 2 + M 2 ^(-M 2 + M 2 ))(4 Mp» - 3 s) 
t D *n*-.D*B* - l.\/j,..U 2 .\/ 2 ' l8J 

For the vector exchange diagrams we obtain the same expression for all the spin states. 
One must do the same calculations for every channel, but, one can see that all the terms 
corresponding to the t-channel have the structure 

(fa + k 3 ) ■ (k 2 + k^tx^u^el 
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Figure 5: S-channel p exchange diagram. 



and those that corresponds to the u-channel 

(h + h) ■ (k 2 + k^tx^iv^,^ 

In Tables El El [KB [III [12 and [Lj of the Appendix we give the D*D* -> channel and 
D*D* — > channel amplitudes projected into isospin and spin. The amplitudes, J/ipJ/tp, 
ujj/ip, 4>J /ip — > K*K* , pp, ujuj, 00, J/i/jJ/if), ujJ/ip, 4>Jip, co>0 for / = and pJ/ip — > K*K* , 
pp, pu, pJ/ip, p<p for 1 = 1, are not in the tables because the interaction is exactly zero 
in our model. Also, the SU(3) amplitudes that involve the K*K*, pp, ujuj, 00 and ujcp 
channels are in the Tables V-X (contact term) and XVIII- XXI (exchange term) in [37] . 
Here we have used the constant g = M p /(2/ 7r ) =4.17 with f n = 93 MeV for light mesons 
and g D = M D ./(2f D ) = 6.9, g Ds = M D ./(2f D .) = 5.47, g Vc = M J/i ,/(2f nc ) = 5.2 with 
f D = 206/ y/2 = 145.66 MeV [36], fn s = 273/^ = 193.04 MeV [36], f Vc = 420/^ MeV 
[46] for the D*, D* and J/tp mesonsj. As we see in these tables the exchange term is equal 
for J = and J = 2 and the exchange of one heavy vector meson is supressed by the factor 
k = Ml/Mjj which is equal to 0.15 in the case of M L = M p and M H = M D *. We can 
observe that the interaction at the D*D* threshold becomes attractive and very strong for 
D*D* -> D*D*, D*D* D*D* and D*D* -> D*D* S amplitudes, which will lead to bound 
states of two vector mesons, basically D*D* and/or D*D*, in the region of interest around 
4000 MeV. 

The formalism that we are using is also allowed for s-channel vector exchange and we 
can have the diagram of Fig. [51 But we found in [34] that this leads to a p-wave interaction 
for equal masses of the vectors, and only to a minor component of s-wave in the case of 
different masses [37], hence, we do not consider these terms. 

-'^We have made the approximation ^gg Va ~ g in the channels that involve one light meson and the J/ip 
meson. 



9 



3 T- matrix 

The results of the amplitudes discussed in the Section 2.2 and 2.3 provide the kernel or 
potential V to be used in the Bethe-Salpeter equation in its on-shell factorized form, 

T = (1 - VG)- X V . (19) 

The potencial V here is a 10 x 10 matrix in I = with the amplitudes obtained from the 
channels 

D*D*(4017), D*D*(4225), K*K*(1783), pp(1551), c^(1565) 

00(2039), J/-0JM6194), cjJ/^(3880), 0J/^(4116), ^0(1802) 

in its elements for each spin J = 0, 1, 2 independently. Also, in/ = l,V A isa6x6 matrix 
whose elements come from the transition potentials between the channels 

D*D*(4017), K*K*(1783), pp(1551), pcj(1558), pJ^(3872), p0(1795). 

In Eq. ffl9l) G is a diagonal matrix where its elements are the two meson loop function Gi 
for each channel i: 

d 4 q 1 1 

(27r) 4 q 2 — m 2 + ie (P — q) 2 — m\ + ie 

which upon using dimensional regularization can be recast as 
„ 1 / m? (Tin — m? + s m 2 , 

= 16^( a + L °V + 2. L '" J 4 

+ JL( Log 3 - '"' + '"■ + 2 ^ + Log 3 ± ~ "i + , (21) 

V s — s + m 2 — mf + 2p v /s — s — m 2 + mf + 2p^/s J J 

where P is the total four-momentum of the two mesons and p is the three-momentum of 
the mesons in the center-of-mass frame: 

J(s- (mi + m 2 ) 2 ) (s - (mi - m 2 ) 2 ) 

p = Vi • (22) 



G < = 1 I ~2 , - T5 ^2 ^2 , ~ » ( 20 ) 



Analogously, using a cut off one obtains 



m a .v 



q dq uoi + u 2 



(,i " ' (2tt) 2 ^ 2 [(P°) 2 - {ui + o; 2 ) 2 + ie] ' (23) 



o 



where q max stands for the cut off in the three-momentum, Ui = (q 2 + m 2 ) 1 / 2 and the 
square of center-of-mass energy (P ) 2 = s. In the complex plane, for a general y/s, the 
loop function in the second Riemann sheet can be written as [47] : 

Gf = Gl(V~s) + Im(p) > (24) 

where G\ refer to the loop function in the second Riemann sheet and G\ is the loop 
function in the first Riemann sheet given by Eqs. (I2T!) and (1231) for each channel i. 
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4 Results 



Once we introduce the potentials obtained from the Tables El El [TUl [HI CE2] and [T3] (in 
the Appendix) as a kernel V in Eq. ( |T9|) . we evaluate the transition matrix T between 
channels and look at \T\ 2 in the real axes and look for poles in the second Riemann 
sheet of the complex plane. If these poles are not very far from the real axis they occur 
in = (M ± iT/2). Since the only meaningful physical quantity is the value of the 
amplitudes for real y/s, only poles not too far from the real axes would be easily identified 
experimentally as a resonance. The amplitude in Eq. ( IT91 close to a pole looks like 

Tij » -^L_ , (25) 

S Sp 

The constants g\ (i = channel), which provide the coupling of the resonance to one par- 
ticular channel are calculated by means of the residues of the amplitudes. 

We have set the parameters of Eq. (1211) as follow: we have fixed the value of \i as 1000 
MeV for all the channels, the subtraction constant for the loops containing SU(3) mesons, 
that we call as a^, is set to —1.65 in order to find the position of the pole / 2 (1275) as 
in [37], for the channels where the two mesons are heavy we put an = —2.07, value that 
is chosen to get the position of the pole found in S = around 3940 MeV, and for the 
rest of the channels, uJ/ip, (fiJ/ip, pJ?p, we also put ai = —1.65. One could do a study 
of the theoretical uncertainties for the positions and widths to account for underlying 
SU(4) symmetry breaking of the model. This study was done in [38] and one finds similar 
results here. In summary, using the freedom of the model to use /„. or fo and My, but 
readjusting the substraction constants to obtain one of the resonances with the same mass 
as the experiment, one finds changes in the masses smaller than 20 MeV and changes in 
the couplings smaller than 8%. 

We find four poles in / = and one pole in I = 1 . The pole positions and the coupling 
constans gi are given in Tables El [31 HI [5] and [6j In / = we have found one pole with 
a mass ~ 3940 MeV for each spin J = 0, 1, 2 and from the couplings we see that they 
couple stronger to D*D*. For J = and J = 1, the width found is very small, T p = 14.8 
and MeV respectively, but for J = 2 the width is comparatively larger T p = 52 MeV. 
Another pole is also found for 1 = and J = 2 with mass M p = 4169 MeV and width 
T p = 132 MeV which couples now stronger to D*D*. In the 1 = 1 sector, we find a pole 
only for J = 2 with mass M p = 3919 and width T p = 148 that couples mostly to D*D*. In 
|21j . the XYZ mesons observed up to now are listed in Table 1. To make the appropiate 
correspondence, let us quote from [21] the properties of the XYZ states discussed, which 
we show in Table [H 

Looking at this table, we see that there are three states with mass around 3940 MeV. 
The identification is subtle. At first one could think that the three states around 3940 MeV 
could correspond to our states found around this mass. However, all the three experimental 
states have C-parity positive, while our state at 3945 MeV has C-parity negative. Hence, 
this state is a prediction and does not correspond to any of the experimental ones. We 
call it Yp (3945) where p stands for prediction. The other two states that appear in Table 
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State 


M (MeV) 


T (MeV) 


jPC 


Decay modes 


Production modes 


Z(3940) 


3929 ± 5 


29 ± 10 


2++ 


DD 


77 


X(3940) 


3942 ± 9 


37 ± 17 


J p+ 


DD* 


e+e" -> J/^X(3940) 


V(3940) 


3943 ± 17 
3914.3±|i 


87 ±34 
33t^ 2 


J p+ 


ujj/ip 


B -»■ KV(3940) 


X(4160) 


4156 ± 29 


139i^ 3 


J p+ 


D*D* 


e+e" -> J/^X(4160) 



Table 1: Properties of the candidate XYZ mesons. 

□ with masses around 3940 MeV are the Z(3940) and the V(3940). The Z(3940) has been 
seen in the spectrum of DD produced in 77 collisions with mass M = (3929 ±5) MeV and 
T = (29 ± 10) MeV and the two-photon production process can only produce DD in a ++ 
or 2 ++ state. From the angular distributions the Belle measurement strongly favors the 2 ++ 
hypothesis, and for this reason we associate this state to our 2 ++ state, which is found with 
mass and width M p = 3922 MeV and T p = 52 MeV. For the V(3940) there are different 
measurements of the mass and width, whereas Babar reports M = (3914. 3^ 3 ; 8 ) MeV and 
T = (33Ti 2 ) MeV, the Belle's values are M = (3943 ± 17) MeV and T = (87 ± 34) MeV. 
Even though the difference in the values of the mass is not large, there is a considerable 
uncertainty in the width. This is the state that we associate to our ++ state, which we 
find with mass M p = 3943 MeV and width T p = 14.8 MeV. The point in favor of this 
assignment is the result of the calculation of T((3943, + [0 ++ ]) — > uJ/iji), that can be done 
straightforward by means the formula: 

T((3943, + [0 ++ ]) - u J/VO = ?l ^f (26) 

with p the momentum of uo in the resonance rest frame. Taking the coupling gy^jj^ = 
(-1429 — i 216) MeV from Table El we obtain T((3943, 0+[0++]) -> uJ/ifi) = 1.52 MeV, 
which is compatible with T(V(3940) — > uJ/ip) > 1 MeV, obtained from the measured 
product of branching fractions B(B -> XV(3940))i3(V(3940) -> uJ/^j) ((7.1 ±3.4) x lO" 5 , 
reported by Belle, and (4.9 ±1.1) x 10~ 5 according to Babar) together with the assumption 
that B(B — ► i^V(3940)) is less than or equal to 1 x 10~ 3 , the typical value for allowed 
B — > K+ charmonium decays. Thus, we find a natural explanation on why this rate is much 
larger than it would be should it correspond to hadronic transitions between charmonium 
states [2"T] . 

The state V(3940) of Tabled] does not decay into uJ/ip [13]. This is a reason not to 
associate it to our (3943, + [0 ++ ]). This state does not fit into our vector - vector scheme, 
suggesting it has a different nature. 
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The state found with J PC = 2++, M p = 4169 MeV and T p = 132 MeV can be clearly 
identified with the X(4160) by the proximity of mass, width and quantum numbers. 

In the next section, we are going to include the DD channel by means a box diagram 
as it was done in [3H [3H1 E7] . 

V~s pole = 3943 + i7A, / G [J^] = 0+[0++] 



D*D* 


D*D* S 


K*K* 


pp ijjuj 


18810-^682 


8426 + zl933 


10 - ill 


-22 + i47 1348 + z234 



00 J/tpJ/lp toJ/l]) <t*J/lp U(f) 


-1000 — il50 417 + i64 -1429 - z216 889 + zl96 —215 - 


zl07 


Table 2: Couplings <?j in units of MeV for I — 0, J = 0. 




V~s pole = 3945 + *0, I G [J PC ] = 0-[l+-] 




D*D* D*D* K*K* pp uu 00 J/*pJ/*P ujj/ij 




18489 -^0.78 8763 + i2 11 - i38 





Table 3: Couplings gi in units of MeV for / = 0, J = 1. 




^ = 3922 + 126, J G [J^]=0+[2++] 




D*D* D*D* K*K* pp ujuj 


21100 -zl802 1633 + i6797 42 + il4 -75 + i37 1558 + il 


821 



00 J/-0J/-0 UjJ/tp 4>J/4> ^(f) 

-904 - il783 1783 + zl97 -2558 - z2289 918 + z2921 91 - z784 



Table 4: Couplings ^ in units of MeV for I = 0, J = 2. 

5 The Dl) decay mode 

5.1 Evaluation of the DZJ-box diagram 

Here we consider the diagrams of Fig. O in order to take into account the DD decay that 
we introduce in the most important channels: D*D* and D*D*. Our starting point is the 
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Ao/e = 41 69 + *66, I G [J PC ] = 0+[2" 



D*D* 




K*K* 


pp U1UJ 


1225 - i490 


18927 - i5524 


-82 + i30 


70 + i20 3 - i2441 



UJ(f) 



1257 + ^2866 2681 + i940 -866 + i2752 -2617-^5151 1012 + ^1522 

Table 5: Couplings gi in units of MeV for / = 0, J = 2 (second pole). 
VVe = 3919 + i74, J G [J^] = l-[2++] 



D*D* 


K*K* pp 


poj pJ/ip 


p<p 


20267 - i4975 


148 - z33 


-1150-^3470 2105 + z5978 


-1067-^2514 



Table 6: Couplings in units of MeV for / = 1, J = 2. 

set of diagrams of Fig. [6l Take the first one, suppose D* = D* + , D* = D*~ and 7r = 7r°, 

used in Section 



the vextex are provided within the same hidden gauge formalism [31 
2, by means of the Lagrangian 

C v ** = -ig(V' i [$,dr$\) . 
For these particles, we have, for the first diagram of Fig. [6j 



|fiH) 4 ^(-^) 4 ^-2^) 



x (k 3 - 2q)"eW(P + h- 2q) a e^(P + h~ 2qfe 
i i 

x 



X 



q 2 — m 2 D + ie (ki — q) 2 — m 2 + ie 

% % 



7r(fei - q) 



(P — q) 2 — m 2 D + ie (k 3 — q) 2 — + ie 

D*(h) D* s 



D* D*„ 



D' 



D{q) 



D 



D(P-q) 

D*(k 2 ) D*(k 4 ) D* 



■*(h - q) k | | K 

D 



D 



D 



d: d* 



(27) 



[28) 



u 



Figure 6: DD-box diagrams for the D*D* and D*D* channels. 
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Using the approximation that all the polarization vectors are spatial, it is possible to write 
the above amplitude as 

d A q 



_ jf {DD) _ 4 f j? A a . a .q ia e MJ®J®Ai) 

11 — 9d J ^tt) 4 ^ j 
1 



q 2 — m 2 D + ie (ki — q) 2 — m 2 + ie 

x 1 1 (29) 

[P — q) 2 — m 2 D + ie (k% — q) 2 — m 2 + ie 

This integral is logarithmically divergent and as in [34J we regularize it with a cutoff in the 
three-momentum of the order of the natural size, for which we take q ma x = 1.2 GeV. The 
results do not change much if one takes a value around this. After performing analytically 
the dq° integral of Eq. f[2"9l . one finds 



y{DD) 



X 




1 



lu — lo d + ie 



- u D + ie P° - 2uj d + ie P° + 2lu d 
(lu 3 + ojjj + 4 ujuj 2 d + Au 2 uj d - k%(j)r\ (30) 



where u = \J q 2 + m 2 , to D = yq 2 + m 2 D , P° = k® + k\. In Eq. f l30|) we can see clearly the 
sources of the imaginary part in the cuts k\{k%) — uj — uju = 0, P° — 2ujo = 0, which give 
rise to the decays D* + — ► tt°D + and D* + D*~ — > D + D~ respectively. After isospin and 
spin projection, we obain 

t (DD,I=0,J=0) _ 45 y(DD) 

t (DD,/=0,J=2) = 18 y(DD) (31) 

and 

t (DD,I=l,J=0) _ g y(DD) 

t (DD,/=l,J=2) = 2 y(DD) ^ ( 32 ) 

where t^- 05 ) is given by Eq. (!30j) after removing the polarization vectors. That is, we have 
found that the DD decay channel contributes only to J = and J = 2. The case J = 1 is 
forbidden, as already found in [34], [38]. This is because the parity of the D*D* system for 
s-wave is positive, which forces the DD system to be in L = 0, 2. Since the D and D have 
no spin, the total angular momentum J is equal to L in this case. Therefore, only the + , 
2 + quantum numbers have this decay channel. 

For the second diagram of Fig. [6], the evaluation of the integral is very similar to that 
of Eq. (129]) . and it becomes the formula of Eq. (1301) except for a factor and changing 



m n — > rriK and to —> u>k, being uok = \/q 2 + rn? Kl that is 
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(DD,I=0,J=0) _ 1 n r~r(DD) 
K K 
(DD,I=0,J=2) _ A r~ADD) 



with 



t^=»> J =» = 4VT^ (33) 



, „ / 1 \ 3 1 / 1 



fADD) 4 8 /*'' maa; J -6 ( 1 \ 



X 



X 



3 157T 2 JO \UJkJ ujd \k® + UJk + UJd , 

1 1 1 



For the third diagram of Fig. [6] we obtain 



ki — ujk ~ ^>d + ie k® — ujk — ^>d + ie P° — 2ujd + ie 

— — (wjc + u*d + 4 ^kujI + 4u 2 K u D - klu D ) (34) 



,(DD,I=0,J=0) 15 sy(DD) 

,(DD,I=Q,J=2) _ 6 ~ {D d) , s 



with 



e» = AjAiLr***'— — 1 — 1 1 



157T 2 Jo UJUJk^D w + k% + UJ + OJd ki + UJk + ^>D 

1 111 



X 



fci — uj k — uj d + ie k® — oj — oj d + ie P° — 2u D + ie P° + 2uj d 
x (2u D (u D + u K ) 2 + u 2 (2uj d + uj k ) +uj{2uj d + uj k ) 2 -2k° 2 uj D ^ . (36) 

As in [38] we use a form factor for an off-shell ir(K) in each vertex, which is 

F(q) = e~£ (37) 

with A = 1.2 GeV [H]. The real part of the potential coming from the DD box is much 
smaller than the real part of the potential coming from the contact plus exchange terms 
as we can see in the representative figure ,Fig. [7J for the D*D* channel. Then, it can 
be neglected and we keep only the imaginary parts ot the DD box diagrams, which are 
plotted in Fig. [HJ In this figure, we see that the most important contribution comes from 
the DD(ttit) box diagram, the DD(KK) and DD(Ktt) contributions being less relevant. 

5.2 Results with 

In Figs. [HI [10] and [UJ we show the results when one introduces the amplitudes obtained in 
Section 2 and the DD box diagrams of Fig. [6] in the Bethe-Salpeter equation Eq. (fl9!) . 
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Real part of the potential for 1=0 and J=0 



u 
Pi 



0F^ 
-500 
-1000 
-1500 
-2000 



DD box 
contact + exch. 



3600 3700 3800 3900 
V7 [MeV] 



4000 



Real part of the potential for 1=0 and J=2 

Oh =^ 



-500 
^- 1000 
£-1500 
-2000 
-2500 



DD box 
contact + exch. 



3600 3700 3800 3900 
V7 [MeV] 



4000 



Figure 7: Real parts of the potential for the D*D* channel. 




Figure 8: Imaginary parts of the DD box diagrams of Fig. [6j 
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£>*£>*(4017) 
4. x 10 13 p : 

3.xl0 13 : : 

H^.xlO 13 : 

l.xlO l3: : 

i 1 . 

3940 3942 3944 3946 3948 3950 
^ [MeV] 

Figure 9: |T| 2 for / = and J = (left), J = 1 (right), in the main channel D*D*. 

The masses and widths in the real axis obtained from these plots are given in Table [7J As 
one can see, the results for the mass and widths of the resonances have not changed too 
much with respect to those of Section 4, the inclusion of the DD diagrams modifies only a 
few MeV the final width. The effect can be seen in Fig. [12] for the resonance in the 1 = 0, 
J = sector. Thus, the main decay channels of our states are the light vector - light vector 
decay modes as K*K* , pp, ujuj, 00, u<p in I = 0, and pu, p<fi for 1 = 1, and the light vector 
- heavy vector decay modes, ujJ/ip for 1 = 0, and pJ/tp in / = 1. 



I G [J PC ] 


Theory 


Experiment 


0+[0++] 

0- [i+-] 

0+[2++] 
0+[2++] 

1- [2++] 


Mass [MeV] 
3943 

3945 
3922 
4157 
3912 


Width [MeV] 
17 



55 

102 

120 


Name 
V(3940) 

'%(3945)" 
Z(3930) 
V(4160) 

'%(3912)" 


Mass [MeV] 
3943 ± 17 
3914.3tii 

3929 ± 5 
4156 ± 29 


Width [MeV] 
87 ±34 

29 ± 10 


jPC 

J p+ 

2++ 
J p + 



Table 7: Comparison of the mass, width and quantum numbers with the experiment. 



6 Conclusions 

We have made a full study of the vector - vector interaction in the C = and S = 
sector using the hidden gauge formalism. The interaction comes from contact terms plus 
vector meson exchange in the t-channel. We have found a strong attraction in the 1 = 0, 
J = 0, 1, 2 and I = 1, J = 2 sectors, enough to bind the vector - vector system. By looking 
for poles in the second Riemann sheet, we have found five resonances, three of which 
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D* D (4017) 



4.x 10 



3. x 10' 
^2.xl0 6 



1.x 10 




3800 3850 3900 3950 4000 4050 



VT [MeV] 



400000 
300000 
R200000 
100000 



P L 
4100 



D S *D S (4225) 




4150 4200 4250 4300 



V7 [MeV] 



Figure 10: |T| 2 for / = and J = 2 in the main channels D*D* (first pole, left) and D* S D* 
(second pole, right). 
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400000 

<n_300000 
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D* D (4017) 





3800 3850 3900 3950 4000 4050 4100 
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Figure 11: |T| 2 for 7=1 and J = 2 in the main channel D*D* 
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Figure 12: |T| 2 for 1 = and J = before and after the inclusion of the DD box diagrams. 
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can be associated with the experimental data: The state (3943, + [0 ++ ]) to the V(3940), 
the (3922,0+ [2++]) to the Z(3930) and the (4157, 0+[2++]) to the A" (4160). There is no 
experimental counterpart for our state (3945, 0~[l + ~]), which is thus a prediction of our 
model. These three states with mass around 3940 MeV are basically composed by D*D*, 
and decay into pairs of light vectors like K*K*, or light vector - heavy vector as ujj/tp. 
Our model predicts another state around 4160 MeV, (4157, + [2 ++ ]), which we identify 
with the X(4160) state in base to the proximity of mass and widths and C-parity. This 
resonance has J PC = 2 ++ and is mostly D*D*. In the 1 = 1 sector, the attraction is 
weak and we find only one resonance in the case of J = 2, (3912, 1 _ [2 ++ ]), the possible 
association of this state to the AT (3940) is unlikely since our state can decay to DD, though 
with small intensity, but this decay is not seen for the X(3940). The width that we obtain 
is also considerably larger than that of the X(3940), T = 120 MeV. 

As we have shown in this work, the region around 3940 MeV is very interesting and 
there could be more resonances not yet seen in this region. The findings of this work 
should motivate the experimentalist to look into this region in the channels that involve 
light vector - light vector or light vector - heavy vector like K*K* and pJ/ip. 
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Contact 



V-exchange 



Total[I[J p ]] 



D*D* -> D*D* 



D*D* -> 



-> if* if* 



D*D* 



D*D* 



PP 



Usui 



D*D* -> 00 

D*D* -> 0J/V> 
-> ui(j) 



D*D* S -> 



D*D S * -> if* if* 



D*D* 



LULU 



d* s d: -> ^ 

Table 8: Amplitude 
Total is obtained at 



6<£ 
2y/2 g D g Ds 

-2^3 gg D 






4 & 
-2V2gg D , 



^V2g Ds g Vc 
o 

o 



gj, (2M^M p 2 +Mj ; (3Mg+M p 2 ))(4Mf„ -3s) 
4A^M 2 M p 2 

g D g Ds (2M% t +2M 2 Dt -3s) 



\ 2 A/ ( ;.» 
gg D (2M 2 Dt +2M 2 K ,~3s) 

V3gg D (2M 2 D +2M 2 -3s) 

gg D (2M 2 ,+2Mg-3s) 



g D g Vc (2M 2 Dt +2M 2 J/ii -3s) 
gg D (2Ml t ,+M 2 /ii +M 2 ~3s) 





gk( M j/^+ M |)( 4M V 3s) 
2M j/4> M 'i 

gg Ds (2M 2 Dt +2M 2 Kf ~3s) 
V2M 2 D , 





gg Ds {2M 2 Di +2M 2 -3s) 
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8 

flI , iflt|c (2Jtf» : +2Af5 / ^-3 S ) 
v^M 2 
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projected in isospin and spin for 1 = and J 
the threshold of D*D*. 
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Contact 



V-exchange 



Total[I[J p ]] 



D*D* -> D*D* 



D*D* -> 



-> if* if* 

£>*/)* -> pp 
D*D* -> cjcj 
-> 00 

-> 0J/V> 

D*D* -> cj0 



£>*£>* -> if* if* 



D*D* S -> J/^J/V 

UjJ/lp 

WW 



9^ 

3y/2g D g Ds 

3ggD 


















g| (2M^M p 2 +Mj ; (3Mg+M p 2 ))(4M2 . -3s) 
9DgD s (2M2»+2M^-3s) 

gg g (2M| )< ,+2M|.,-3s) 
2M^ 





3k( M j/^+ M |)( 4M or 3s) 

9S b 3 (2M2 j+ 2M|,»-3 S ) 













-46.1 
-23.7 g D g Ds 
-1.3 gg D 










-10-3 

2Agg Ds 





Table 9: Amplitude 
Total is obtained at 



projected in isospin and spin for / = and J = l.The approximate 
the threshold of D*D*. 
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Contact 



V-exchange 



Total[I[J' 



D*D* -> D*D* 



D*D* -> 



D*D* 
D*D* 
D*D* 
D*D* -» 
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-> pp 

— > CiJCJ 

+ <i>J/rl> 



D*D* -> cj0 
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"2 & 
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2 Mi 





9 D9t)c (2M| ) »+2Mj /;/j -35) 
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g<yi5 s (2Af| ) ,+2Af2. H ,-3s) 





sflD3 (2M2 J+ 2M|-3 S ) 
5Ds5 „ c (2M2 J+ 2Jkf2 /Vj _3 s ) 



gg g3 (2M^+M^+M|-3s) 







Table 10: Amplitude projected in isospin and spin for / = and J 
Total is obtained at the threshold of D*D*. 



-58.1 ^ 
-29.4 ^ Ds 
5-4 ^ D 

10.1 

-5.8 gg D 


-7.2 ^ D ^ c 
9.5 ^ 



-18.3^. 




-7-3gg Ds 
-4.6 ^ Ds ^ c 


8-5 gg Ds 


2. The approximate 
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Contact 


V-exchange 


~ Total[I[J p ]] 


D*D* - 


> D*D* 


2 02, 


»», (2M2M p 2 +M5 /i/) (-M,;+M p 2 ))(4M| ) »-3 S ) 


0.6 g 2 D 


4M^M 2 M 2 


D*D* -> 


K*K* 


~99d 


gg D (2M^+2M^-3s) 

S 




D*D* 


-> pp 











D*D* 


— > pcj 


-2y/2gg D 


g g D (2M 2 *+M 2 +M 2 -3 S ) 
V2M 2 Df 




D*D* - 




^99D 


gg D (2M 2 Dt +M 2 J/jl; +M 2 p -3s) 


-3.5 gg D 


D*D* 














Table 11: Amplitude projected in isospin and spin for I — 1 and J = 0. The approximate 
Total is obtained at the threshold of D*D*. 





Contact 


V-exchange 


~ Tota/[/[J p ]] 


D*D* -> 

-> pp 
_D*D* -> puj 
D*D* -> pJ/> 
D*D* -> p0 


9ggD 


S 2 , (2M2M 2 +M 2 /V) (-M2+M p 2 ))(4M2 , -3s) 
4M 2 MlMl 

J ftp w P 

gg D (2M%, +M%, (4M 2 ? +2M|,» -3s)+2M 2 , (2M 2 -3s)) 


1-602, 
-2.5 gg D 

-1-3 ggD 





2 

3<7<?D 


9 9c (2M 2 < ,+2M 2 -3s) 


72 






2v / 2Af 2 » 







Table 12: Amplitude projected in isospin and spin for 1 = 1 and J = 1. The approximate 
Total is obtained at the threshold of D*D*. 
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Contact 


V-exchange 


^Total[I[J p ]] 


D*D* - 


> D*D* 


-9h 




-2.452, 




D*D* -» 


K*K* 


99D 
2 


SSD (2Af|,»+2M-2.»-3s) 

S 


2-7 ^ D 


D*D* 


-»• pp 











D*D* 


— > pcu 


V2gg D 


gg D (2M*,+M*+M*-38) 
V2M 2 D , 




D*D* - 




-ZggD 


ggn(2M% sf +M^+M^-3 S ) 


-9.5 gg D 


D*D* 














Table 13: Amplitude projected in isospin and spin for 7 = 1 and J = 2. The approximate 
Total is obtained at the threshold of D*D*. 



27 



